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Abstract 

Analyses of betatron coupling can be broadly divided into two categories: the 
matrix approach that decouples the single-turn matrix to reveal the normal modes 
and the hamiltonian approach that evaluates the coupling in terms of the action 
of resonances in perturbation theory. The latter is often regarded as being less 
exact but good for physical insight. The common opinion is that the correction of 
the two closest sum and difference resonances to the working point is sufficient to 
reduce the off-axis terms in the 4x4 single-turn matrix, but this is only partially 
true. The reason for this is explained, and a method is developed that sums to 
infinity all coupling resonances and, in this way, obtains results equivalent to the 
matrix approach. The two approaches is discussed with reference to the dynamic 
aperture. Finally, the extension of the summation method to resonances of all 
orders is outlined and the relative importance of a single resonance compared to all 
resonances of a given order is analytically described as a function of the working 
point. 

1 Introduction 

The example of linear betatron coupling will be used in the first instance to demonstrate 
the summation of the influences of all the resonances in a given family into a single driving 
term. It will then be shown that is in fact a general result that can be applied to all linear 
and non-linear resonances. 

Analyses of betatron coupling can be broadly divided into two categories: the matrix 
approach ||, Q, || that decouples the single-turn matrix to reveal the normal modes 
and the hamiltonian approach ||, [|7| that evaluates the coupling in terms of the action 
of resonances using a perturbation method. The latter is often regarded as being less 
exact but good for physical insight. The general belief is that the correction of the two 
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closest sum and difference resonances to the working point should be sufficient to reduce 
the off-axis terms in the 4x4 single-turn matrix, but in most cases this is not successful. 



1.1 Matrix method for coupling compensation 

The 4x4 single-turn matrix in the presence of skew quadrupoles and/or solenoids is of 
the form 



(where M,n,m,N G 3? 2x2 ). Coupling compensation is achieved by setting the two 2x2 
matrices n and m to zero. Due to symplecticity and periodicity of T only four free 
parameters (that is the strengths of four compensator units) are required. However, this 
compensation is only valid at the origin of T. 

A transformation can also be applied to the matrix T that decouples the linear motion, 
so making it possible to describe the beam in the whole machine with the well-known 
Courant and Snyder parametrisation in the transformed coordinates. 

1.2 Classical hamiltonian method for coupling compensation 

This method is based on the expansion in a Fourier series of the coupling perturbation 
term in the Hamiltonian. The standard procedure is to assume that the low-frequency 
components dominate the motion and that only the nearest sum and difference resonances 
therefore require compensation (single resonance compensation). 
The essential difference between this and the matrix approach is: 

• The matrix method is exact while the hamiltonian method is approximate; 

• A coupling compensation made by the matrix method is only valid at one point in 
the ring whereas the hamiltonian method gives a global correction; 

• The matrix method leaves finite excitations in all resonances, including those closest 
to the working point, whereas the hamiltonian method leaves finite excitations only in 
the far resonances. 

The reason for the two last points is that the matrix method includes all resonances 
automatically and combines them in such a way that the matrix is uncoupled at one 
point, while the hamiltonian method sets only the closest sum and difference resonances 
to zero. If the far resonances have little effect, then the two methods are virtually equiv- 
alent. This is however an uncommon situation. 

The logical implication is that by finding a way to sum all the resonances, the classical 
hamiltonian method can be made to reproduce the results found with the matrix method. 
Once this is done, the natural questions are which of the two methods is the better for 
operation, and if the principle can be extended to higher orders. 




(1) 
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The aims of this paper are the following: 

• To outline a summed resonance compensation procedure (taking into account both 
the low and high-frequency part of the perturbative Hamiltonian) for the case of linear 
coupling and to extend this result to the non-linear case (Section 0); 

• To analytically compare (Section |3|) the single and summed resonance theories point- 
ing out some general results that can be obtained using the analytical expression of the 
generalized driving term; 

• To numerically compare the single and the summed resonance compensations for the 
linear coupling (the latter is shown to be equivalent to the matrix compensation) using a 
4D coupled Henon map f§ (Section |); 

2 Multiple resonance compensation for linear cou- 
pling 

The starting point for this analysis is taken from ||, with the initial assumptions that: 

• the perturbative Hamiltonian is calculated at 9 = 0. Since the origin is an arbitrary 
choice, this is not a restriction; 

• solenoid fields are absent. The presence of solenoid fields does not change the argu- 
ment but by omitting them the resulting equations become more transparent. They will 
be added at the end. 

The linear coupling compensation using the notation of || requires (without any ap- 
proximation for removing the high frequency part): 

' +oo 

^ioio-p = sum resonance 

p=— oo 

< ' (2) 

+oo 

ftjooi-p = difference resonance. 

k p=-oo 

2.1 Detailed derivation for difference resonances 

Consider first the treatment of the difference resonances and express ([|) explicitly as 

+ 00 OO „27T 

E - ^oo = E / M0)e-W* -Qv- 9 (3) 

p=— oo P=— oo 

where 

A W) = " ^ m H9), (4) 
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9 = s/R is the coordinate along the ring, Q x ,y are the horizontal and vertical tunes, /j, x y 
are the horizontal and vertical phase advances, /3 x>y are the horizontal and vertical beta 
functions and R is the radius of the ring. 

Suppose k(6) is different from zero and constant in j short intervals (i.e. the regions 
occupied by the sources of coupling and by possible correctors) [Oi, Qi + AQi] in which A{6) 
can be considered approximately constant (thin lens approximation) q 

C 'oo = E(A^oo)*- ( 6 ) 

i=l 

Now define A~ = Q x — Q y and compute (AC~)$, the contribution to from the i — th 
sub-element, 

(AO= E / A(0) e -^ A -#d0, (7) 



p=—oo 

which, when integrated gives, 



(ACJ, ~ Aft) [ e ~ l[ ~ Pli + ] ~ e ~ t[ ~ Pl \ ■ ( 8 ) 

p=— oo 

The summation can be redefined making use of the shift [A - ] (the closest integer to A - ) 
so that p = [A~] + k where k is an integer. In the limit (A - — [A - ]) A6i -C 1, the previous 
expression becomes^] 

+oo 

(AC-), * A^e-^' ~ [ADft J2 (A - _ rl-n _ + M) ~ e m 

k=-oo ^ t " 

= A(O i )e- i ( A ~ " [ A "D e « £ - {sin(^) - S in[fe(^ + Aft)] + 

,±4 (A -[A ])-* 

+ z [cos[k(9i + A9i)} - cos(k9i)]}. (9) 

To sum these series first we rewrite them in a more suitable form (x = d{ or x = Q\ + AQ,[) 
valid when < x < 2tt |J: 



+oo 



\- sin(fcx) _ 7rsin[(A - [A ])(% - x)] _ 

(A- - [A-]) - * sin((A- - [A-])*) ' 1 ) 



1 In a real machine A(0) will vary slowly or, at least, it will be possible to cut the elements into short 
enough pieces that A{&) can be considered as constant over all sub-elements to any desired degree of 
accuracy. 

2 In the expression (^) the term e~ l( - A ~' A ^ Ae is expanded to the zero order while the terms e lfeAe 
with k ~ (A - — [A - ]) are not expanded. This assumption, supported "a posteriori" by the accuracy of 
the final result, is based on the fact that the contribution of the higher order terms to the sum of the 
series is negligible. 
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cos(£;x) 7rcos[(A —[A — x)] 

k Z^ oo (A--[A-])-k = sin((A- - [A-])tt) ' ( } 

The application of ( |T0D and ( ]TT] ) to (||) then gives 

(AC"), = A(6 t )e-^ ~ [A1) ^ sin((A _ z [A - ])7r) (( sin [( A ' ~ [ A D(* ~ g > ~ A ^ + 

- sin[(A- - [A-])(tt - 0*)] + i[cos[(A- - [A~))(it - 9, - A0 4 )] + 

- oos[(A--[A-])(7r-e < )]]} = 

-i[(A~ - [A-])(7r-0< -)] 

■ e LV 1 JA 2 ;J . (12) 



After expressing A(^) explicitly, (A(7 co )i becomes 

fc, (A--[A-])A^ 
4i2sin((A- - [A-])tt) 



fAC-V ~ ^ (A--[A-])Ag, / 3(e)3(c) J M) - toffl - (A- - [A-])tt) 



(13) 



Equation (|T3|) can be summed directly for all the elementary elements in the ring to 
give the coupling coefficient for the combined influence of all the linear difference 
resonances. 

2.2 Extension to sum resonances and solenoids 

For the sum resonances the procedure is unchanged and the formal result is the same, 
but with the following substitutions: 

- Hy^i) ► fJt x (6i) + Hy{Bi) (14) 

A- — > A+ = Q x + Q y . (15) 



In presence of a uniform solenoidal field the summed resonance driving term can again be 
written in the form 

+ 00 „27T i 

Ct= / ^(^ _l(A ~ p)9 d9 (16) 
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where now 
A±{9) = 

(k — 0, S - 



AixR 
R 



P X (9 i )Py(9 i )RS 



— I 



1 r 1 

Px Py 



2Bp 

(AC+), and (ACJ,: 



Px Py . 

Bq). The same procedure as before yelds the following expressions for 



(ACS* 



,S t sin((A ± -[A ± ])A^/2) 
' 2 sin((A± - [A±])tt) 



l l 



z[^)±// 2/ (^)-(A ± -[A ± ])7r] 



(18) 



2.3 Extention to the nonlinear case 



It is clear from Section ^]2| that the procedure for summing the resonances is, in fact, 
independent of the detailed form of the term A(9) and that, with the general form of 
A(9) from ]16| , the method can be extended to the nonlinear case. 

The driving term of a given resonance (ni,n 2 ) of order N = n\ + n 2 for the single and 
summed resonance theories are (respectively) 



1 r 

nun2 ' p ~ 7r(2R) N / 2 \n 1 \\\n 2 \\J Px Py 1 ' 



and 



n - [AD i o|m|/2o|n 2 |/2 i[n 1/la: +n 2Aty -(A-[A])7r]r,j/) 

Un ^°°- sm[n(A-[A])]n(2R) N ^\n 1 \\\n 2 \\J Px Py 

(20) 



whereQ: 



1 j 2 2\Bp\ [[ } 2 [ dx( N -M- 2 )dy(M+V dx^-M-^dyM* 
d^B s a x . . d^B s m n 2 >, 

■(Kl^r- - Mtt) - l ^rau^(y - tt)1 ( 21 ) 



3 In the following formulas the partial derivatives are evaluated in x = y = 
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for |ri 2 1 even, N > 3 and 1 < | 1 < {N — 2); 

. ( it! 2 jn^ dV-VBy fr N -*B x 

^ ' 2\Bp\^ 2 * , &r(AH"al- 2 )0y(l«al+ 1 ) dx^-^-^dy^V 

. 0l x \ a y\ 

for \n 2 \ odd, N > 3 and 1 < |n 2 | < (JV - 2). 

It is customary to use the symbols C for the coupling driving terms and K for the 
higher-order non linear driving terms. There is also often a factor of 2 between the two 
definitions (K — C/2 for a given resonance). Since this report concentrates on coupling 
the "C" -styled definition has been extended to all the cases. 



3 General results 

This section is dedicated to pointing out some of the general consequences of the summed 
resonance theory. 

3.1 Analytic comparison of the influence of single and summed 
resonances 

It is interesting to compare the contribution to the coupling excitation from all resonances 
to that of the closest single resonance. The single resonance driving term reads 

C p = A(6)e- l ( A -P^d6 (23) 
Jo 

where now p is the closest integer to A. 

Using the thin lens approximation the i — th contribution to C p can be written 

-A6i 



(AC p )i ~ A(6i) / l 'e-i^-P^M 
hi 



A{6l) \A±-p\ e ~* {A ~ P)6% [C ° S[(A " p)Ml] Sin[( A " p)Mt]] ~ 

k ' A0 < (^)i ni|/2 (/^)i^ l/ V[ ni ^^ +n3 ^-t A - 1 '^]. (24) 



7r(2 J R) JV / 2 |ni|!|n 2 |! 



Fig. ([I]) shows the ratio between ||(AC*)j|| and ||(ACp)j|| varying the distance from the 
resonance ([A 1 * 1 ] = 1). 
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Figure 1: Ratio r between ||AC^|| and ||AC^|| versus A ± . 



The formulae ( ]T5| ) and (23) give the same result for the modulus of the driving terms 
when exactly on resonance. A difference between (|13|) and (|24| ) increases (approxima- 
tively) quadratically as one moves away from A integer. Agreement between the two 
formalism can be therefore expected only if the working point is close enough to the res- 
onance to be compensated []. 

However, this is not usually the case if the aim is the full compensation of the linear cou- 
pling, that is both the sum and difference resonances. In most pratical cases, the working 
point is close to the difference resonance and relatively distant from the sum resonance. 



3.2 Closed-orbit distortion from a dipole kick 

The equation ( p~3| ) can be applied to the resonance family 

Qz=P (25) 
where z = x,y. This leads to the expressions for the closed-orbit distortion due to a dipole 



kick and links the single resonance theory ||, |T6[ to the integrated theory of Courant 
and Snyder [17] . 

4 Note that the phase terms are different even when exactly on resonance. 
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In this case 

hoo „2tt 



Coo = V / A^e^^-^)^ 



(26) 



where now [16 



= (27) 

(AB/Bp is the dipole error). 

For a localized error of length Al (thin lens approximation): 

(AC )• = - Qz ~ [Qz] ^***e i M0)-(Q'-[Q'W = 

{ °° >l " 23/2 J Ri/2 S m( 7 r(Q a -[QJ)) V/3 * 



1 g z - [gj AlAB 

23/2^1/2 sin ( 7r g ;s ) 5p 



rp-jinM-Qz*). (28) 



Comparing equation (|28|) to the closed-orbit distortion (at the origin) due to a kick oc- 
curring at a given position 8 

(z)e=a = V^M^ cos(/*,(0) - vrQ z ) (29) 

Zsm[7TLJ z ) np 

shows that the normalized orbit distortion differs from Re(AC 00 )j by only a constant, 

1 Re(AC 00 ) J (30) 



%j e=0 2-vm-v\Q z -[Q z \) 

3.3 Betatron amplitude modulation 

Applying the same procedure to the resonance family 

Qz = 2p. (31) 

one gets the modulation of the betatron function due to a small gradient error occurring 
at a given position 6. 
In this case 

+ 00 „27T 

Coo== J2 / 4(0)eW*-p)0<10 (32) 

p=-oo ^ 

with 

^ W = J_ ftWe ^W|^. (33) 
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This gives 



(AC )■ = -- Qz - [Qz] B(9) — ^e^M-2(Q z -iQzM = 

[lXUcoh " 2sin(27r(Q z - [Q z ])f z[U) Bp dx 

= 4^V#1^W#^( 2 ^W- 2 ^). (34) 
2 sm(27r<y z ) Bp ox 

The last expression coincides with the modulation of the beta function (at the origin) 

< AA >« = l^kj^wt coM) ~ Q - T)) (35) 

and shows that the normalized modulation differs from Re(ACoo)j by only a constant, 

'A&\ 1 



3.4 Comments 



A;„-(ft-[Qj)75 Re(ACJ - (36) 



The summed resonance driving terms lead naturally to a new definition of bandwidth. 
However, pratical differences maybe not so clearly observed. When the working point is 
close to a resonance, the bandwidth is important but the single and summed theories do 
not differ much. When the working point is far from the resonance, the bandwidth is 
unimportant and any widening may go unnoticed. 

There is also the more academic point that the magnitude of the summed resonance 
driving term is dependent on the azimuthal position in the machine. The closed orbit is 
a very good example of this. If there is a closed bump in the orbit, then the Re(C , i i0iO o) is 
zero outside the bump and finite inside. The implication is that the beam is responding to 
standing waves from each member of the resonance family. The summed response is the 
sum of these standing waves. When inside a bandwidth the growth rates are related to 
the standing wave amplitude and varying according to the position around the machine. 



4 The coupled Henon map 

In this section, the summed resonance approach is shown to be equivalent to the matrix 
approach and both are compared to the single-resonance compensation by performing 
a numerical analysis on the so-called Henon map || : a hyper-simplified lattice model [] 
whose phase-space trajectories show some of the expected characteristic of a realistic lat- 
tice map (nonlinearities, regions of regular and stochastic motion etc.) In this application 
the linear coupling is generated and corrected by 1+4 thin skew quadrupoles 0. 

5 A linear lattice model containing only one sextupolar kick. 

6 Lattices with only solenoids or with both types of coupling elements give the same kind of results. 
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The global compensation of the coupling resonances (at #=0) is achieved if 

( 5 

^[Rfi(AC-)+iIm(AC-)] = 

i=i 

(37) 

5 

^[Re(AC+) + iIm(AC+)] = 0. 

The compensation for both the sum and the difference resonance is obtained solving the 
system of 4-equations (for the 4 unknowns ki) given by (|37|) . 



k (m 2 ) 


Matrix 


Summed 


Single 


k\ (source) 


0.5 


0.5 


0.5 


k 2 


-0.051 


-0.050 


0.559 


h 


0.034 


0.033 


0.554 


ki 


-0.319 


-0.313 


0.476 


h 


-0.275 


-0.275 


0.117 



Table 1: Compensator strengths (fe-s) i n presence of the coupling source k\ using the 
single-turn matrix compensation and the summed and single resonance compensations. 

Table 1. shows a comparison between the strengths of the 4 correctors (/C2-5) when com- 
pensating the single-turn matrix, the two infinite families of sum and difference resonances 
(for the same 8 = 0) and the closest sum and difference resonances to the working point. 
The single-turn matrix compensation has been performed by means of the MAD program 
fllOf while the (single and multiple) resonance compensation has been obtained making 
use of the AGILE program Jll]] in which the formula (|i~3"D has been implemented. 
The single-turn matrix (in 9 = 0) in presence of the coupling source ki = 0.5 m~ 2 (no 
compensation) has non zero off-axis 2x2 sub-matrices given by 

0.49 0.49 \ 



/ M 


n ^ 




\ m 


Nj 





0.02 0.02 



-0.23 
-0.01 



-0.24 
-0.01 



(38) 



while the residual values of n and m after the single resonance compensation {C + 
0) are given by 

/ -0.11 -4.67 \ 

0.02 0.15 

0.03 -4.57 
\ 0.03 -0.09 / 



c- 



(39) 
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The off-axis terms are in fact larger after the compensation then before. This is explained 
by the influence of the far resonances that can not be neglected, to get a satisfactory 
coupling compensation in this case. 

The same conclusion may be drawn by looking at the driving terms of the closest sum 
and difference resonances before the correction 

|C + | = 0.0628 |C~| = 0.0628 (40) 

and afterwards 

\C + \ = 0.0207 \C~\ = 0.1018. (41) 

Last two equations show that the sets of the uncompensated sum and difference reso- 
nances have a "weight" comparable (larger in the case of the difference resonance) to the 
ones compensated. 



The quantitative difference between the two approaches can be better investigated by 
means of a tracking analysis. 

In the following, the results from stability and footprint diagrams as well as the calculation 
of the dynamic aperture for the compensated Henon map are shown. 



4.1 Stability and footprint diagrams 

A stability diagram and the related frequency diagram can be obtained by the follow- 
ing procedure: for each initial condition inside a given grid in the physical plane (x, y) 
(Px — Py — 0), the symplectic map representing the lattice is iterated over a certain 
number of turns. If the orbit is still stable after the last turn, the nonlinear tunes can be 
calculated using one of the methods described in . In the stability diagram the stable 
initial conditions are plotted whereas in the frequency diagram are represented the corre- 
sponding tunes. The insertion in the frequency (footprint) diagram of the straight lines 
representing the resonant conditions up to a certain order makes it possible to visualize 
the excited resonances close to the reference orbit. 



Figs. (5) - (10) show the stability and frequency diagrams for the uncoupled Henon 
map, after the single resonance compensation and after the summed resonance one. 
The comparison points out that the summed resonance compensation allows a more ef- 
ficient restoration of the uncoupled optics. It is significant the analysis of the degree of 
excitation relative to the resonances (3,-6), (1,-4) and (2,-5) for the two different compen- 
sation approaches. 



Using the perturbative tools of normal forms |13j one can calculate the value of the 
first resonant coefficient (leading term) in the interpolating Hamiltonian for the consid- 
ered resonances. The leading term can be considered as a "measure" of the resonance 
excitation. It can be shown [14| that in absence of coupling the leading term of the res- 
onances (3,-6) and (1,-4) is different from zero (first order excitation) whereas the one of 
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the resonance (2,-5) is zero (second order excitation). The strength of the coupling (that 
is, in the considered case, the strength of the residual coupling after the compensations) 
is proportional to the growth of the leading term of the first order non-excited resonances 
and to the decrease of the leading term of the other ones. 

The resonance degree of the excitation varying the compensation approach can be better 
visualized ploting the network of the resonances and their widths inside the stability do- 
main. The analysis of Figs. (2)-(4) confirms that the SR method is characterized by a 
residual coupling considerably stronger than the one left by the MR compensation. 
The same conclusion can be drawn the following topological argument. A trace of the 
presence of linear coupling in a nonlinear lattice is the spliting of the resonant channels 
in correspondence of the crossing points (multiple resonance condition in the tune space). 
This phenomenon is evident only in the case of the MR compensation (see the central 
part of Fig. (4)). 
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Figure 2: Network of resonances of the uncoupled Henon map. 
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Figure 3: Network of resonances after the summed resonance compensation. 



4.2 Dynamic aperture calculations 



The dynamic aperture as a function of the number of turns N can be defined |L5[] as the 
first amplitude where particle loss occurs, averaged over the phase space. Particle are 
started along a grid in the physical plane (x, y): 



x = r cos 9 



y = r sin 9 



(42) 



and initial momenta p x and p y are set to zero. 

Let r(9, N) be the last stable initial condition along 9 before the first loss (at a turn 
number lower than N). The dynamic aperture is defined as 



D 



i9,N)]Um(29)d9 



(43) 



An approximated formula for the error associated to the discretization both over the radial 
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Figure 4: Network of resonances after SR compensation. 



and the angular coordinate can be obtained replacing the dynamic aperture definition with 
a simple average over 9. Using a Gaussian sum in quadrature the associated error reads 



AD 



(Ar) 



+ 



dr 
89 



(A6>) 5 



(44) 



where Ar and A9 are the step sizes in r and 9 respectively. 

In Tab. 2 the values (with the associated errors) of the dynamic aperture are quoted for 
the three studied optics for short (iV=5000) and medium (iV=20000) term tracking. 

The difference between the summed and the single resonance compensations is notice- 
able: the compensation of the all families relative to the coupling resonances allows an 
improvement close to 10% respect to the case in which the high frequency part of the 
perturbative hamiltonian is neglected. 

It can also be pointed out that the summed compensation seems to slightly improve (at 
the limit of sensitivity due to errors) the dynamic aperture respect to the uncoupled case. 
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D (m) 


Uncoupled 


Summed 


Single 


iV=5000 


0.0406 


0.0412 


0.0372 


iV=20000 


0.0405 


0.041 


0.037 



Table 2: Dynamic aperture values relative to the uncoupled Henon map and after the 
summed and single resonance compensations. The associated error (according to the 
formula 0) is about 2% for iV=5000 and about 4% for iV=20000. 

5 The nonlinear case 

Dealing with high-order resonances with a view to optimizating stability is not so straight- 
forward as in the linear case: the number of resonances that can be excited both by a 
given multipole and by the set of correctors meant for compensating a given resonance 
becomes higer and higher; moreover the resonance compensation is only one of the tools 
that has to be used to get a succesfully optics optimization. 

For these reasons we have not here attempted a general comparison between the 
summed and the single resonance compensations using tracking analysis. We intend to 
return to this question in the future. 

We note however that a certain number of attempts to compensate one particular sex- 
tupolar resonance for the Henon map show that the two compensations are not far only 
if the working point is close enough to the considered resonance. The summed resonance 
compensation is to in general better in the case of the compensation of several resonances 
at the same time. 
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6 Conclusions 



A general method has been derived for the summation of all the resonances within a given 
family both for the linear and for the non-linear cases. The fact that this summation is 
valid and gives a meaningful result is confirmed by its application to the known closed- 
orbit distortion equation, the betatron modulation equation and the decoupling of the 
linear transfer matrix for a ring. The application of the summed-resonance driving term 
to the coupling raises the question of the relative merits of the different types of coupling 
compensation that are now possible. This problem has been investigated with the help 
of the Henon map. The results indicate that use of the summed-resonance compensation 
(equivalent to the matrix approach) yields a larger dynamics aperture. 
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Figures caption: 



• Figure 5: Stability domain of the uncoupled Henon map. 



• Figure 6: Footprint diagram of the uncoupled Henon map. 



• Figure 7: Stability domain after the summed resonance compensation. 



• Figure 8: Footprint diagram after the summed resonance compensation. 



• Figure 9: Stability domain after the single resonance compensation. 



• Figure 10: Footprint diagram after the single resonance compensation. 
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